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ABSTRACT 


The  purpose  of  this  paper  is  to  introduce  three  discrete  distributions  named  Poisson  exponential  distribution, 
Poisson  size  biased  exponential  distribution  and  size  biased  Poisson  exponential  distribution.  These  distributions  apply  to 
biological  data  sets,  traffic  datasets  and  thunderstorm  datasets.  These  distributions  are  introduced  with  some  of  its  basic 
properties  including  moments,  coefficient  of  skewness  and  kurtosis  are  discussed.  The  method  of  moments  and  maximum 
likelihood  estimation  of  the  parameters  of  discrete  PED,  PSBED  and  SBPED  are  investigated.  It  is  found  that  the 
reciprocal  of  MOM  and  MLE  estimator  is  unbiased  for  the  proposed  distributions.  Applications  of  the  three  models  to 
different  discrete  data  sets  are  compared  with  Poisson  distribution,  size  biased  Poisson  distribution,  size  biased  generalized 
Poisson  distribution  size  biased  geometric  distribution  and  size  biased  Poisson  lindley  distribution  to  test  their  goodness  of 
fit  and  the  fit  shows  that  the  proposed  distributions  can  be  an  important  tool  for  modelling  biological,  traffic  and  other 
discrete  data  sets. 


KEYWORDS:  Goodness  of  Fit,  Estimation  of  Parameters,  Exponential  Distribution,  PED,  Poisson  Distribution, 
PSBED,  SBPD  Moments,  Size-Biased  Exponential  Distribution,  Size-Biased  Poisson  Distribution 

INTRODUCTION 


Poisson  distribution  was  derived  by  French  mathematician  Simeon-Denis  Poisson  in  1830with  probability  density 
function  (pdf) 

p-Aj* 

P(X)  = -^L,  * = 0,1,2 (1.1) 

Where  X > 0 is  the  shape  parameter 

The  exponential  distribution  is  the  distribution  of  continuous  random  variables  with  probability  density  function 

(pdf) 

/(*)  = Xe~Xx,  x > 0 (1. 2) 

Where  X > 0 is  the  scale  parameter 

Weighted  distribution  was  first  introduced  by  Fisher  (1934)  and  was  later  formalized  by  Rao  (1965).  When  the 
probability  of  selection  of  an  event  is  not  equal  weights  are  assigned  to  the  distribution  to  remove  bias. 

The  pdf  of  weighted  distribution  is 


fw  (*:60 


w[x)f0Cx- 9) 

£[»W] 


(1.3) 
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When  w(x)  — xr  The  weighted  distribution  are  called  moment  distributions.  For  r= 1 and  2 moment  distribution 
is  named  as  size-biased  and  area-biased  distributions  simultaneously. 

When  w (x)  =x,  the  weighted  distribution  is  called  size  biased/length  biased  distribution 

f(x;  6)  = xfo(x:d)  (1.4) 

v- 1 

E(x)  =n\ 

Weighted  distribution  has  a number  of  applications  in  environmental  sciences,  including  ecology  and  forestry  and 
fields  of  real  life  including  medicine,  ecology  and  reliability. 

The  pdf  of  the  size  biased  Poisson  distribution  is 

afr/Q  = 1 Sr  (L5) 

X=l,  2,  3...  A > 0 

The  pdf  of  the  size  biased  exponential  distribution  is 

f{x ; 0)  = x92e~xe  (1.6) 

X=l,  2,  3...  0 > 0 

The  mixture  of  two  or  more  probability  distributions  is  known  as  mixture  distribution.  The  probability  distribution 
is  called  a mixture  distribution  if  its  pdf  is  written  in  the  form 

P(X  = x)  = f f(x/A)gA(X)dA  (1.7) 

Shakila  Bashir  and  Mujahid  Rasul  (2016)  introduced  a discrete  distribution  named  Poisson  are  a biased  Lindley 
distribution  and  its  application  to  biological  data.  Rama  Shankar,  Hagos  Fasshaye,  Abrehe  Yemane  (2015)  applied  the 
SBPLD  to  model  thunderstorm  data.  They  proved  that  SBPLD  gives  much  closer  fit  than  SBPD  and  gives  a better 
alternative  to  the  SBPD  and  can  be  recommended  for  modelling  thunderstorm  events. 

Shankar  and  Mishra  (2015)  discussed  a size-biased  two  parameter  Poisson  Lindley  distribution(SBTPPLD)  by 
size  biasing  the  two  parameter  Poisson  Lindley  of  Shankar  and  Mishra  (20 14) which  has  been  obtained  by  mixing  Poisson 
distribution  with  two  parameter  Lindley  distribution  introduced  by  Shankar  and  Mishra  (2013).  Bhatti  (2015)  derived  the 
new  generalized  Poisson  Lindley  distribution. 

Porinita  Dutta  and  Muhindra  Borah  (2014)  further  investigated  some  properties  and  applications  of  SBPLD, 
such  as  moments,  cummulants,  co-efficient  of  variation,  harmonic  mean,  survival  function,  estimation  of  parameters  by  the 
method  of  moments  etc.  Khurshid  Ahmed  Mir  and  Munir  Ahmed  (2009)  introduced  size-biased  discrete  distribution, 
generalized  size  biased  discrete  distribution  and  their  generalization.  Ghitani  et  al  (2008)  obtained  the  size-biased  Poisson 
Lindley  distribution  and  its  properties.  Shaban  (1981)  derived  the  Poisson-inverse  Gaussian  distribution,  by  mixing  the 
Poisson  with  the  inverse  Gaussian  distribution.  Bulmer  (1974)  discussed  the  Poisson  lognormal  distribution  by  mixing  the 
Poisson  distribution  with  lognormal  distribution.  Sankaran  (1970)  derived  poison-Lindley  distribution  by  mixing  Poisson 
with  Lindley  distribution.  He  discussed  estimation  of  parameters  and  applied  data  on  proposed  distribution.  He  proved  that 
Poisson  gives  a better  fit  as  compared  to  some  other  distributions. 
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In  this  paper,  we  consider  the  Poisson-exponential  distribution  PED,  PSBED  and  SBPED  which  is  obtained  by 
integration  and  discretization  of  the  mixture  of  Poisson  and  exponential  distribution,  Poisson  and  Size-biased  exponential 
distribution  and  Size-biased  Poisson  and  Size-biased  exponential  distribution. 

1.  POISSON  EXPONENTIAL  DISTRIBUTION 

The  Poisson  exponential  distribution  (PED)  arises  from  the  Poisson  distribution  with  pdf  (1.1) 

f(x/A)  = x = 0,1,2, ... 

When  its  parameter  A follows  the  exponential  distribution  with  pdf  (1.2) 
g(A)  = 0e~eA  A = 0,1,2, ... 

So 

P(X  = x)  = f f(x/X)gx(A)dA 
P(X=x)  = /*— — 6e~BX  dA 

*'0  x\ 

After  simplification  the  pdf  of  the  PED  is 

P(X=x)  = ^4^’  6 > 0,x  = 0,1,2 (2.1) 

PROPERTIES  OF  POISSON  EXPONENTIAL  DISTRIBUTION 


The  rth  factorial  moment  of  the  PED  is 


• _ rYr 

^ (r)  — gr 

For  r=l,  2,  3&  4 the  first  four  factorial  moments  of  PED  are 

' _ i ' _ 2 ' _6  ' _24 

^ (1)  - e’  ^ (2)  “ 02  ’ M (3)  - i?  ’ ^ (4)  - ^ 


Therefore  the  first  four  raw  moments  of  PED  are 


■ _ 1 _ CF  ' _ 2 + 0 - 

Mi  g X,  g.  2 g2’M3 


82+6B+6  • _ 03+1402+360+24 

03  ’ M 4 _ 04 


The  mean  moments  of  PED  are 
l+e 


Mz  = 


02 


(B+2)(0+l) 


P4 


0 +100  +180+9 


The  co-efficient  of  skewness  and  kurtosis  of  PED  are 

(0  + 2)2 


Pl  = - 


0+1 


(2.2) 


(2.3) 
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of  e. 


Pl  = 


e3+ioe2+i8e+9 


02+20+l 

The  moment  ratios  for  PED  show  that  for  different  values  of  0>O,  pi  increases  and 
For  0=0,  p2=9  and  as  0— P2  is  increasing. 

The  moment  generating  function  of  the  PED  is 


Mx(t)  = 


e 


1+0-e1 

The  Cummulant  generating  function  and  the  four  cummulants  of  PED  are 

(0+2)(9+l)  ,,  03+702+120+6 


K1  = i K2=^,K3  = - ' , K4  = 

The  probability  generating  function  of  the  PED  is 

, x 0 

§(t)  = 

The  Fisher  information  matrix  of  the  PED  is 
F = 


92(i+9) 

Some  more  properties  of  PED  are 


2 

(1  — G 


e 2 


(2.4) 


(2.5) 


(2.6) 


(2.7) 


It  follows  that  for  all  values  of  6 > 0 the  poisson-exponential  distribution  is  over-dispersed  unless  for  large  values 


/(x+l;0)  _ {0|(l+0)*+2} 

/fe0)  ~ {0|(l+0)*+1} 

/(x+l;0)  _ 1 

/(x;0)  _ 1+0 


Figure  1.1:  Graph  of  p & a2  of  PED 


Figure  1.2:  Graph  of  p & pt  of  PED 


Impact  Factor  (JCC):  3.9876 


NAAS  Rating  3.45 


Family  of  Poisson  Distribution  and  its  Application 


5 


Figure  1.3:  Graph  of  p & p2  of  PED 


• Figure  (1.1)  shows  the  trend  of  mean  and  variance  of  PED.  As  0 increases  the  mean  and  variance  of  PED 
decreases,  so  p and  c2  have  the  same  trend. 

• Figure  (1.2)  shows  the  trend  of  mean  and  coefficient  of  skewness  of  PED.  As  0 increases  the  mean  decreases  and 
coefficient  of  skewness  increases. 

• Figure  (1.3)  shows  the  trend  of  mean  and  coefficient  of  kurtosis  of  PED.  As  0 increases  mean  of  PED  decreases 
and  coefficient  of  kurtosis  for  PED  increases. 

• For  the  PED  model  (f^  , fi2)  — > (4,  9)  as  0->O 

• For  the  PED  model  the  moment  ratios  increase  as  0 increases. 

ESTIMATION  OF  PARAMETERS 

In  this  section  the  parameter  of  the  PED  is  estimated  by  MOM,  MLE  and  properties  of  the  estimates  are  derived. 
Method  of  Moments  (MOM) 

The  method  of  moment  estimate  (9 of  the  parameter  0 is  given  by 


X 


Maximum  Likelihood  Estimates  (MLE) 

The  maximum  likelihood  estimate  9 of  the  parameter  0 is  given  by 
9 = - 

X 

Theorem  1 

If  X follows  PED  with  pdf  P(x)  = —ztt  the  MOM  and  MLE  estimator  = is  unbiased  for  - and  as  n -> 

F v 2 (0  + l)*+1  9 9 

co  var  Q-)  — > 0 

Proof 

The  pdf  of  the  Poisson  - exponential  distribution  is 

~ ce+i)x+1 

0 = j 9^es  | = X(A) 
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Applying  expectation  on  (A) 


= E(X) 


After  simplifications 


= E{X) 


Since  the  mean  of  the  PED  is  - 

e 


Ew=l 

i i 

Which  proves  that  ~ is  unbiased  for  - i.e  the  reciprocal  of  the  estimator  is  unbiased. 
0 0 


var  (j-J  = i mr(X) 


After  simplification 

var  (j)  = i var(X ) (f?) 


l+e 


Variance  of  the  PED  is  var  (X)  =— p putting  this  expression  in  equation  (B)  on  RHS 


var 


© 


i+e 

n02 


Hence  lim,,^  var  ( ~ ) = 0 


Which  proves  that  the  variance  of  the  PED  is  minimum  as  n ->  oo 

2.  POISSON  SIZE  - BIASED  EXPONENTIAL  DISTRIBUTION 


The  Poisson  size-biased  exponential  distribution  arises  from  the  Poisson  distribution  with  pdf  (1.1) 
f{x/X)  = x = 0,1,2, ... 

When  the  parameter  A follows  the  size-biased  exponential  distribution  (SBED)  with  pdf  (1.6) 
f{pr,  6>)  = x02e~xe 
So 

PiX  = x)  = f f(x/A)gA(A)  dA 
P(X  = x)  = foK(^f-)A02e-ASdA 
After  simplifications  the  pdf  of  PSBED  is 

p(*  = *)  = a£x(x  + 1> 

0>O,  x=0,  1,2,3 


(2.8) 


(2.9) 


(3.1) 
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PROPOERTIES  OF  POISSON  SIZE  - BIASED  EXPONENTIAL  DISTRIBUTION 

The  rth  factorial  moment  of  PSBED  is 


fr+2 


(r)  gr 


(3.2) 


For  r=l,  2,  3 & 4 the  first  four  factorial  moments  of  PSBED  are 


. _ 2 . _ 6 , _ 24 

A*  (i)  " e ’ 11  (2)  “ iz  ’ A*  (3)  - p-  ’ A*  (4) 


120 

04 


Therefore  the  first  four  raw  moments  of  PSBED  are 


. _ 2 r _ 6+20  __ 

A1 1 — g ’ A4  2 ~ ^2  ’ A4  3 - 


24+180+202 


03 


= ' 


203+4202 +1440+120 
04 


The  mean  moments  of  PSBED  are 
20+2 


M2 


Fs  = 


^4  = 


02 


202+60+28 
~03 

203+2602+480— 72 
04 


The  co-efficient  of  skewness  and  kurtosis  of  PSBED  are 

404+2403+14802+3360+784 


Pi 

P2 


803+2402+240+8 

03+1302+240— 36 
202+40+2 


(3.3) 

(3.4) 


The  moment  ratios  for  PSBED  show  that  for  different  values  of  0>O,  pi  increases  and  for  0=0,  p2=98  and  as  0-+° 
p2  is  increasing 


The  moment  generating  function  of  PSBED  is 
02 

02+20-(20+2)e£+ezt+l 

The  probability  generating  function  of  PSBED  is 
ft)  = e2 

^ 2 02  + 20+t2— 2t(0  + l)+l 

The  Fisher  information  matrix  of  PSBED  is 
F = 


02(i+0) 

Some  more  properties  of  PSBED  are 


(3.5) 


(3.6) 


(3.7) 


A i = °- 


02 


It  follows  that  for  all  values  of  0 > 0 the  poisson  size  biased  exponential  distribution  is  over-dispersed  unless  for 
large  values  of  9. 
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/(x+l;0)  _ 92(2x+l)/(9+l)x+3 
f^x-, 9 ) ~~  02(x+l)/(0  + l)*+2 

gx+i-,9)  = (2+yx) 

fix-, 9 ) (1+Vz)(e  + 1) 

This  shows  that  is  a decreasing  function  in  x,  that  is  f(x;0)  is  logarithmically  concave  function. 

Therefore  PSBED  is  unimodal. 


Figure  2.1:  Graph  of  p & a2  of  PSBED 


Figure  2.2:  Graph  of  p & fa  of  PSBED 


Figure  2.3:  Graph  of  p & p2  of  PSBED 

• Figure  (2.1)  shows  the  trend  of  mean  and  variance  of  PSBED.  As  0 increases  the  mean  and  variance  of  PSBED 
decreases,  so  p and  o2  have  the  same  trend. 

• Figure  (2.2)  shows  the  trend  of  mean  and  coefficient  of  skewness  of  PSBED.  As  0 increases  the  mean  decreases 
and  coefficient  of  skewness  also  decreases,  therefore,  mean  and  coefficient  of  skewness  has  the  same  trend. 

• Figure  (2.3)  shows  the  trend  of  mean  and  coefficient  of  kurtosis  of  PSBED.  As  0 increases  mean  of  PSBED 
decreases  and  coefficient  of  kurtosis  for  PSBED  increases. 

• For  the  PSBED  model  (ft  ,ft)  ^(98,  -18)  as  0->O 
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ESTIMATION  OF  PARAMETERS 

In  this  section  the  parameter  of  PSBED  is  estimated  by  MOM,  MLE  and  properties  of  the  estimates  are  derived. 
Method  of  Moments  (MOM) 

The  method  of  moment  estimate  6 of  the  parameter  0 is  given  by 

S = 2- 

x 

Maximum  Likelihood  Estimates  (MLE) 

The  maximum  likelihood  estimate  8 of  the  parameter  0 is  given  by 

e = 2- 

x 

Theorem  2 


If  X follows  PSBED  with  pdf  P(x)  = 


02(x+l)  1 1 

— — the  MOM  and  MLE  estimator  - is  unbiased  for  - and  as  n 

(0+i)*+2  e 0 


S) 


co  var  ~ 


Proof 


The  pdf  of  Poisson  size  biased  exponential  distribution  is 

02(x+i) 


P(x ) = 


(0+i)*+2 


x 2 . 1 X , 

9 =j  gives  - = -(A) 


Applying  expectation  on  (A) 

E 6) = s (!) 

After  simplifications 


E ~ 


(!)  = ;«« 


Since  mean  of  PSBED  is  - Therefore,  putting  this  value  in  equation  above 


E\s)  = b 

1 1 

Which  proves  that  ~ is  unbiased  for  - i.e  the  reciprocal  of  the  estimator  is  unbiased. 
0 0 


(3.8) 


var  Q-)  = var  (0 
After  simplification 

var  (I)  = (S) 
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Variance  of  PSBED  is  var  (X)=  2 , putting  this  expression  in  equation  (B  ) on  RHS 

m (i+0) 
var  ~ = — -- 
\e)  2nd2 

Hence  lim)W00  var  Qj  = 0 

Which  proves  that  variance  of  PSBED  is  minimum  as  n ->  oo 


(3.9) 


3.  SIZE  - BIASED  POISSON  EXPONENTIAL  DISTRIBUTION 

The  size-biased  Poisson  exponential  distribution  (SBPED)  arises  by  giving  weights,  W(x)  =x  to  Poisson 
exponential  distribution,  Using  the  pdf  of  the  weighted  distribution  from  (1.4) 

f(x;  0)  = xfo<<x’B^  with  E(x)  = ' 
i 1 1 

The  pdf  of  size-biased  Poisson  exponential  distribution  is 

v 


0 > 0,  X=l,  2,  3...00 


Above  is  the  pdf  of  SBPED  which  can  be  derived  by  integrating  mixture  of  size  biased  Poisson  and  size  biased 
exponential  distribution  from  (1.5)  and  (1.6)  also 

The  pdf  of  size-biased  Poisson  distribution  (SBPD)  is 


g{x/X)  = 


e~^Ax~1 


(x— 1)! 


When  its  parameter  A follows  a size-biased  exponential  (SBED)  model  with  pdf 

h(A,6)  = 62Ae~dl 

then 


P{x)  — f0'  g(x/A)h(A,  0)dA 

= C ‘L-?—  o2Ae~0XdA 

J0  (x— 1)! 

After  simplifications  the  pdf  of  size-biased  Poisson  exponential  distribution  is 

nr  'i  x02 
PM  = (TTfl)VU 

9 > 0,X=1,  2,  3, °° 


(4.1) 


Therefore,  when  Poisson  exponential  distribution  was  given  weights  its  pdf  is  the  same  when  we  made  a mixture 
of  SBPD  and  SBED. 
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PROPERTIS  OF  SIZE  - BIASED  POISSON  EXPONENTIAL  DISTRIBUTION 


The  first  four  raw  moments  of  size-biased  Poisson  exponential  distribution  are 

' _ 0+2  ' _ 02+60+60  ' _ 03+1402+360+24  . _ 04+3O03+15O02+24O0+12O 

Mi  - — > M2  - p ’ M 3 - p ’ M 4 — p 

The  mean  moments  of  PSBED  are 
26+2 


F-2  = 


92 


Ms 


(U4  = 


202+360+4 

e3 

24 +480 +860 2— 580 3 
04 


The  co-efficient  of  skewness  and  kurtosis  of  SBPED  are 

404+ 1440 3 +13 120 2 +2 880 +16 


Pl  = 


803+2402+960+8 
12 +240 +43 02— 290 3 


20z+40+2 

The  moment  ratios  for  SBPED  show  that  for  different  values  of  0>O,  pi  increases  and  for  0=0,  P2=6 

The  moment  generating  function  of  SBPED  is 

..  , , , 92et 

~ 02+20-(20+2)e£+ez£+l 

The  probability  generating  function  of  SBPED  is 
, , e^t 

^ 2 02  + 20+t2— 2t(0  + l)+l 

The  Fisher  information  matrix  of  SBPED  is 
F = 


(4.2) 

(4.3) 


02(i+0) 

Some  more  properties  of  SBPED  are 


(4.4) 


(4.5) 


(4.6) 


/r  - er2  = 1 + 


0 02 


8 4 

03  04 


It  follows  that  for  all  values  of  9 > 0 the  size -biased  Poisson  exponential  distribution  is  over-dispersed  unless  for 
large  values  of  0. 

f(x+l-,9 ) 02(x+l) 


/(x;0)  (1 +0)*+1+1 

/O+i;0) 


The  result  above  shows  that  ^ x+1,e l is  a decreasing  function  in  x,  therefore  PSBED  is  unimodal. 

/(*;0)  S 
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Figure  3.1:  Graph  of  p & a2  of  SBPED 


Figure  3.2:  Graph  of  p & P-1  of  SBPED 


Figure  3.3:  Graph  of  p & P-2  of  SBPED 

• Figure  (3.1)  shows  the  trend  of  mean  and  variance  of  SBPED.  As  0 increases  the  mean  and  variance  of  SBPED 
decreases,  so  p and  a2  have  the  same  trend. 

• Figure  (3.2)  shows  the  trend  of  mean  and  coefficient  of  skewness  of  SBPED.  As  0 increases  the  mean  decreases 
and  coefficient  of  skewness  increases. 

• Figure  (3.3)  shows  the  trend  of  mean  and  coefficient  of  kurtosis  of  SBPED.  As  0 increases  mean  of  SBPED 
decreases  and  coefficient  of  kurtosis  for  SBPED  also  decreases. 

• For  the  SBPED  model  (ft  , ft)  ->(2,  6)  as  0->O 

ESTIMATION  OF  PARAMETERS 

In  this  section  the  parameter  of  SBPED  is  estimated  by  MOM,  MLE  and  properties  of  the  estimates  are  derived. 
Method  of  Moments  (MOM) 

The  method  of  moment  estimate  6 of  the  parameter  0 is  given  by 

8 = — 

x-l 

Maximum  Likelihood  Estimates  (MLE) 
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The  maximum  likelihood  estimate  6 of  the  parameter  0 is  given  by 

8 = — 
x-l 


Theorem  3 

X 0 ^ 11 

If  X follows  SBPED  with  pdf  P(x ) = — — — — the  MOM  and  MLE  estimator  ~ is  unbiased  for  - and  as  n —> 

r J (e+i)x+1  e e 

oo  var  Q)  ->  0 

Proof 


The  pdf  of  size  biased  Poisson  exponential  distribution  is 

nr  -v  _ x92 

' ' ~ (0+1  )*+1 


~ 2 . 1 x— 1 , 

0 = _ gwes_  = _(A) 


Applying  expectation  on  (A) 

£©  = £(Jr) 


After  simplifications 


The  mean  of  SBPED  is  [A  = — — Putting  in  the  equation  above 

0 

After  simplifications 


Which  proves  that  =■  is  unbiased  for  \ i.e  the  reciprocal  of  the  estimator  is  unbiased. 
0 0 

VCLV  (f)  = VaV  (~2 ~) 

var  (i)  = mr(x)  ( B ) 


20  + 2 

Variance  of  SBPE  distribution  is  var(X)=  , putting  this  expression  in  equation  (B)  on  RHS 


var 


© = 


(i+g) 

2nd2 


Hence  var  (~  ) = 0 


Which  proves  that  variance  of  SBPED  is  minimum  as  Tl  — * oo 


(4.6) 


(4.7) 
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Table  1 


Measures 

PD 

PED 

PSBED 

SBPED 

A, 

1/ 

2/ 

0 + 2 

I1 

'6 

/0 

0 

T 

1 + 0 

20  + 2 

20+2 

CN 

02 

02 

iPi 

lyVx 

b + 2)2 

|04  + 603  + 3702  + 840  + 196 

40 4 + 14403  + 131202  + 2880  + 16 

J 0 + 1 

J 203  + 602  + 60  + 2 

803  + 2402  + 960  +8 

P2 

3 + 1/i 

03  + 1O02  + 180  + 9 

03  + 1302  + 240  - 36 

24  + 480  + 8602  _ 5803 

02  + 20  + l 

2(0  + l)2 

402  + 80  +4 

The  table  above  shows  the  mean  (|i),  variance  (o2),  coefficient  of  skewness  (i//?i)  and  coefficient  of  kurtosis 
(/?2)  for  the  Poisson  distribution  (PD),  Poisson  exponential  distribution  (PED),  poisson  size  biased  exponential 
distribution  (PSBED)  and  size  biased  Poisson  exponential  distribution  (SBPED). 

APPLICATIONS 

To  illustrate  the  applications  and  to  justify  the  suitability  of  the  three  distributions  i.e:  Poisson  exponential 
distribution,  Poisson  size  biased  exponential  distribution  and  size  biased  Poisson  exponential  distribution  in  practical 
applications,  we  are  fitting  these  distributions  to  some  datasets. 

GOODNESS  OF  FIT  FOR  POISSON  EXPONENTIAL  DISTRIBUTION  (PED)  AND  POISSON  SIZE  - 
BIASED  EXPONENTIAL  DISTRIBUTION  (PSBED) 

In  this  section  we  have  applied  Poisson  exponential  distribution  (PED)  and  Poisson  size  biased  exponential 
distribution  (PSBED)  to  biological  data  and  traffic  datasets  using  the  MOM  and  MLE  estimates  and  compared  with 
Poisson  distribution  (PD). 

A:  Shankar,  et  al  (2015)  gave  data  on  Observed  and  expected  number  of  Hemocytometer  yeast  cell  counts  per 
square  observed  by  ‘student’  1907. 

Form  table-2  it  can  be  seen  that  PED  and  PSBED  give  much  closer  fit  than  the  PD  to  the  yeast  cell  data  set, 
therefore  PED  and  PSBED  are  suitable  for  modelling  data  in  ecology. 

B:  Shankar,  et  al  (2015)  gave  data  on  Observed  and  expected  number  of  red  mites  on  apple  leaves. 

From  table-3  it  can  be  seen  that  PED  and  PSBED  gives  a better  fit  than  PD,  therefore  PED  and  PSBED  are  can  be 
considered  as  an  important  tool  for  modelling  datasets  in  ecology. 

C:  Gerlough,  D.L  (1955)  studied  the  Analysis  and  comparison  of  accident  data.  He  gave  data  on  the  observed  and 
expected  number  of  car  accidents. 

From  table-4  it  can  be  seen  that  PED  and  PSBED  seems  to  be  a good  replacement  of  PD  to  model  traffic  data  sets 
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Table  2:  Goodness  of  Fit  for  PD,  PED  and  PSBED  on  Number  of  Hemocytometer  Yeast  Cell  Counts 


Number  of  Cells  Per  Square 

Observed  Frequency 

Expected  Frequency 

PD 

PED 

PSBED 

0 

128 

118.1 

128.0838 

123.6145 

1 

37 

54.3 

40.35407 

46.22106 

2 

183 

12.5 

12.71395 

12.96199 

3 

3 

1.9 

4.005654 

3.231105 

4 

1 

0.2 

1.262021 

0.755095 

5+ 

0 

0.0 

0.397612 

0.169404 

Total 

187 

Estimation  of  parameters 

0 = 0.459893 

0=2.174386 

0=4.348838 

x2 

9.903 

2.966108 

3.38411 

d.f 

1 

2 

1 

p-value 

0.0016 

0.226944 

0.0658 

Table  3:  Goodness  of  Fit  for  PD,  PED  and  PSBED  on  Number  of  Red  Mites  on  Apple  Leaves 


Number  of  Mites  Per  Leaf 

Observed  Frequency 

Ex 

pected  Frequency 

PD 

PED 

PSBED 

0 

38 

25.3 

37.2101 

32.22676 

1 

17 

29.1 

19.90271 

23.52404 

2 

10 

16.7 

10.64544 

12.8786 

3 

9 

6.4 

5.693964 

6.267186 

4 

3 

1.8 

3.045552 

2.859222 

5 

2 

0.4 

1.628986 

1.25226 

6 

1 

0.2 

0.871302 

0.53322 

7+ 

0 

0.1 

0.466036 

0.222415 

Total 

80 

80 

79.46 

79.76 

Estimation  of  parameters 

0=1.15 

0=0.8696 

1.739 

X2 

18.275 

2.398823 

4.82913 

d.f 

2 

3 

2 

p-value 

0.0001 

0.4939 

0.0894 

Table  3:  Goodness  of  Fit  for  PD,  PED  and  PSBED  on  Number  of  Car  Accidents 


Number  of  Car  Accidents 

X 

Number  of  Roads  Observed 

fn 

Expected  Frequency 

PD 

PED 

PSBED 

0 

18 

14.5 

21.09371 

18.33412 

1 

14 

16.4 

11.20605 

13.26296 

2 

7 

9.3 

5.953225 

7.19585 

3 

4 

3.5 

3.162656 

3.470334 

4 

1 

1.0 

1.680163 

1.569032 

5 

0 

0.2 

0.892588 

0.681026 

6 

0 

0.04 

0.474188 

0.287383 

7 

1 

0.007 

0.251913 

0.118796 

8 

0 

0.001 

0.133829 

0.04834 

Total 

45 

45 

44.8483 

44.9678 

Estimation  of  parameters 

0 =1.1333 

0 =0.88235 

0 =1.76471 

x2 

3.84 

1.2922 

1.3881 

0.05733 

d.f 

1 

2 

2 

P-value 

0.25561 

0.49955 

0.9718 

GOODNESS  OF  FIT  FOR  SIZE  - BIASED  POISSON  EXPONENTIAL  DISTRIBUTION 

In  this  section  we  have  tried  to  fit  size-biased  Poisson  distribution  (SBPD),  size-biased  geometric  distribution 
(SBGD)  and  size-based  Poisson  exponential  distribution  (SBPED)  to  different  datasets  like  thunderstorm  data,  number  of 
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workers  N,  has  I accidents  and  number  of  migrants  aged  15  years  and  above  (survey  1978  data)  using  the  MOM  and  MLE 
estimates. 

A:  Shankar,  et  al  (2015)  gave  data  on  Observed  and  expected  frequency  of  thunderstorm  events  containing  X 
thunderstorm  at  Cape  Kennedy  for  August. 

Table-4  shows  that  SBPED  gives  much  closer  fit  than  the  SBPD  and  SBPLD  to  the  thunderstorm  dataset. 
Thus  SBPED  can  be  considered  as  an  important  tool  for  model  thunderstorm  data. 

B:  Mir,  K A.,  Ahmed,  M (2009)  gave  data  on  Observed  and  expected  frequencies  of  the  number  of  workers  N, 
having  I accidents. 

From  table-5  it  can  be  seen  that  SBPED  gives  much  closer  fit  to  model  count  datasets  than  SBPD  and  SBGPD, 
therefore  SBPED  provides  a better  alternative  to  SBPD  and  SBGPD  to  model  count  data  sets. 

C:  Priti.,  Singh,  B.  P.  (2015)  gave  data  on  Observed  and  expected  number  of  households  with  at  least  one  male 
migrant  according  to  the  number  of  male  migrants  aged  15-years  and  above  (survey  1978  data)  in  the  semi  urban 
household. 

Table-6  shows  that  SBPED  gives  a better  fit  than  SBGD  and  can  be  considered  as  an  important  tool  to  model 
migration  data  sets. 


Table  4:  Goodness  of  Fit  for  SBPD,  SBPLD  and  SBPEDs  on  Thunderstorm  Event 


X 

Observed  Frequency 

fo 

Expected  Frequency 
fe 

SBPD 

SBPLD 

SBPED 

1 

201 

194.2 

199.5 

199.663 

2 

60 

68.3 

59.9 

59.68365 

3 

10 

12.0 

13.3 

13.38056 

4 

3 

1.4 

2.6 

2.666496 

5 

2 

0.1 

0.7 

0.498171 

Total 

276 

276.0 

276.0 

275.8919 

Estimation  of 
parameters 

6 = 0.351449 

9 = 6.37477 

0 = 5.69071 

r 

1.414 

0.165 

0.15495 

d.f 

1 

1 

1 

p-value 

0.2344 

0.6846 

0.6939 

Table  5:  Goodness  of  Fit  for  SBPD,  SBGPD  and  SBPED  on  Number  of  Workers  Having  Accidents 


Number  of  Accidents 

Number  of  Workers 

fo 

Expected  Frequency 

fe 

SBPD 

SBGPD 

SBPED 

1 

2039 

2034.27 

2039.83 

2039.078 

2 

312 

319.48 

309.76 

311.2761 

3 

35 

33.45 

36.38 

35.63846 

4 

3 

2.63 

3.66 

3.626934 

5 

1 

0.17 

0.37 

0.346044 

Total 

2390 

2390.0 

2390.00 

2389.965 

r 

0.772 

0.069 

0.0111245 

d.f 

1 

1 

1 

p-value 

0.379599 

0.7928 

0.916093 
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Table  6:  Goodness  of  Fit  for  SBGD  and  SBPED  on  Number  of  Households  with  at  Least  One  Male  Migrant 


Number  of  Migrants 

Observed  Frequency 

fo 

Expected 

Frequency 

r« 

SBGD 

SBPED 

1 

95 

86.54 

86.52835 

2 

19 

31.32 

31.32305 

3 

10 

8.50415 

4 

2 

11.14 

2.05232 

5 

3 

0.464334 

Total 

129 

128.8722 

Estimation  of  parameters 

e = 0.819 

6 = 4.5249 

f 

7.01 

3.1621 

d.f 

1 

1 

p-value 

0.0082 

0.075366 

CONCLUSIONS 

The  PED,  PSBED  and  SBPED  are  discrete  distributions  that  are  obtained  by  a mixture  of  Poisson  and 
Exponential  distribution  (PED),  a mixture  of  Poisson  and  Size  biased  exponential  distribution  (PSBED)  and  by  applying 
weights  on  PED.  Some  important  properties  of  the  above  distributions  are  derived.  From  Figure  (1),  Figure  (2)  and  Figure 
(3)  it  can  be  seen  that  these  distributions  are  positively  skewed.  The  estimation  of  parameters  is  done  by  the  method  of 
moments  (MOM)  and  maximum  likelihood  estimate  (MLE)  and  it  is  proved  that  the  reciprocal  of  MOM  and  MLE 
estimator  is  unbiased  for  PED,  PSBED  and  SBPED.  The  discussion  on  estimation  and  applications  of  PED,  PSBED  and 
SBPED  demonstrates  that  the  above  distributions  give  a better  approach  to  the  biological  data  sets,  traffic  datasets  and  is  a 
good  alternative  to  PD,  SBPD,  SBGD,  SBGPD,  SBPLD. 
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